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Abstract 

This paper presents the first numerical solution to the non-linear evolution equation for 
diffractive dissociation processes in deep inelastic scattering. It is shown that the solution 
depends on one scaling variable r = Q"^ /Q^'^{x,xq), where Qf (x,xo) is the saturation scale 
for the diffraction processes. The dependence of the saturation scale Q^{x,xq) on both x 
and xq is investigated, (Yq = ln(l/a;o) is a minimal rapidity gap for the diffraction process). 
The X - dependence of turns out to be the same as of the saturation scale in the total 
inclusive DIS cross section. In our calculations Q^{x,xq) reveals only mild dependence on 
xq. The scaling is shown to hold for x <^ xq but is violated at a; ~ xq. 
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1 Introduction 



Diffractive inclusive production in deep inelastic scattering (DIS) at high energy has become an 
area of particular interest of experts since it provides a deeper insight into dynamics of QCD in 
the kinematic region where the density of partons is expected to be high (see Ref. 0] and reference 
therein) . 

Inclusive diffraction in DIS offers an opportunity to probe the transition region between "soft" 
and "hard" interactions giving natural estimates for the value of the shadowing corrections in 
DIS, namely AF2 = F2 — F^'^^^p = ^ as was firstly shown in Ref. [Q] on the basis of the 
AGK cutting rules A more detail approach started with the Kovchegov-McLerranj^ formula 
which expresses the ratio of the diffraction cross section {(Tdiff) to the total cross section (atot) in 
DIS initiated by the quark-antiquark pair produced in 7* ^ g + g decay of the virtual photon. 
This formula reads 



R 



ad^ff _ J dHj dzj d''r^P^''{z,rr,Q')N\r^,x;b) 



atot 2jd^JdzJ d^r^ P^*(^,r^;Q2) N{r^,x;b) ' 

where N{rj_, x; b) is the imaginary part of the elastic dipole-target amplitude for dipole of the size 
r± scattered at fixed Bjorken x = Q^/W^ {Q"^ is the photon virtuality and W is its energy in the 
target rest frame) and at fixed impact b. P^* {z, r_|_; Q"^) is the probability to find a quark-antiquark 
pair with size inside the virtual photon ||^, |^: 
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^ Z} { {z^ + (1 - zf)a^Kl{ar^) + AQ^ z\l - zfKl{ar^) }, 
/ 



where 0? = z{l — z)Q'^ + m^. The functions \&t,l stand for transverse and longitudinal polarized 
photon wave functions. Eq. (|1 . 1| ) is important since it provides a relation between the dipole- 
target elastic amplitude and the cross section of the diffraction dissociation. A non-linear evolution 
equation was derived for the former [|, |1^, [11], |12|, [1^, This equation has been studied both 
analytically 0, |T§ and numerically [g [T^, 0. 

The formula (|1 . 1|) fails to describe correctly the experimental data on the diffraction production. 
Moreover, inclusion of an extra gluon emission in the initial virtual photon wave function is still 
insufficient to reproduce the data [0, 20, 21, 22 1. Nevertheless, Eq. (|1 . 1|) can be viewed rather as 



initial condition to a more complicated equation. 

The non-linear equation for the diffraction dissociation processes can be written for the amplitude 
A^^ which has the following meaning p3 . 



We introduce the cross section for diffraction production with the rapidity gap larger than given 
Yo = ln(l/xo): 



adiffix,Xo,Q^) = J d^r±J dz P^* {z,r^;Q^) a^^^i^{r±,x,Xo) 



;i.3) 



^Fd is the diffractive structure function introduced in Ref. 



1 



and 

dipole 



(^SLir±,x,Xo) = d%N^{r^,x,Xo;b). (1.4) 



The function A^^ is the amphtude of the diffraction production induced by the dipole with size 
r_L with rapidity gap larger than given (lo). Note that the minimal rapidity gap Yq can be 
kinematically related to the maximal diffractively produced mass: xq = {Q"^ + M'^)/W'^. 

The non-linear evolution equation for A^^ was derived in Ref. |^ and recently rederived in Ref. 
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7r"= Jyo 

[ d2x2^^[2 iV^(xo2, y, Yo; b - ^xi^) + iV^(xo2, y, Y^; h - ^xi2)iV^(xi2, y, Y^; h - ^xos) 

-4 Ar^(xo2, 2/, l"o; b - ^xi2)iV(xi2, y\h - ^xo2) + 2 Ar(xo2, y\h - ^xi2)Ar(xi2, y\h - ^xo2)] • 

(1.5) 

The evolution (pTSp is a subject to initial conditions sX x = Xq. 

N'^ir^, xo, xo; b) = N\r^, x^- h) . (1.6) 

Namely, at the energy equal to the energy gap diffraction is purely given by the elastic scattering 



as it was stated in Eq. (1.1 



Since at high energies color dipoles are correct degrees of freedom we can write the unitarity 
constraint : 

2A^ = A^^ + F, (1.7) 
where the function F denotes contributions of all the inelastic processes. An important obser- 



vation is that F satisfies the same equation as A^ [11, 12 1 but with shifted initial conditions 
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Another interesting quantity to study is the cross section of diffractive dissociation process with 
a fixed gap or equivalently to a fixed mass: 

3fJ = - dN^/dYo . (1.9) 



The function 3? was introduced in Ref. |^. The authors of this paper proposed a model in 
which 3? was shown to possess a maximum when varying Yq at fixed Y. Physically this maximum 
means that at given Y there is a preferable mass for the production. Below we will argue that the 
appearance of the maximum is related to the scaling phenomena to be displayed by the function 
AT^. 
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The present paper is entirely devoted to the numerical solution of the equation (|1.5|) . Various 
properties of the solutions A^^ are investigated while our final goal computation of the diffraction 
cross section will be published separately ||2^. In the next Section (2) the solution of the equation 
( |1.5| ) is presented. Section 3 deals with the determination of the diffractive saturation scale. 
Scaling phenomena is discussed in Section 4. We conclude in the last Section (5). 



2 Solution of the non-linear equation 

In this section we report on the numerical solution of the equation ( |1.5|) . The method of iterations 
proposed in Ref. [0 is applied. The constant value for the strong coupling constant as = 0.25 
is always used. The solutions are computed for 4 x 10~^ < a^o < 10~^ and within the kinematic 
region 10~^ < a; < xq and distances up to a few fermi. 

The function A^^ is formally a function of four variables: the energy gap Xq, the Bjorken variable 
X, the transverse distance r^, and the impact parameter h. The 6-dependence is parametric only 
because the evolution kernel does not depend on h. In order to simplify the problem we will 
proceed similarly to the treatment of the 6-dependence of the function In that paper 

we assumed the function to preserve the very same 6-dependence as introduced in the initial 
conditions: 

N{r^,x]h) = (1 - e-"("'''^)^('')), (2.10) 
with the function n being related to the "6 = 0" solution A^(r_|_, x): 

K(x,rx) = - ln(l - N{r^,x)). (2.11) 



N{r±,x) represents a solution of the very same non-linear equation (see Refs. jlll, 0) but with 
no dependence on the third variable. The initial conditions for the function N{rj_,x) are taken 
at 6 = 0. For the case of the proton target ||l6l the anzatz in the form (|2.10|) was shown to be 



a quite good approximation of the exact 6-dependence of the solution to the non-linear equation 
for N{r±, x; b). In Ref. [T^ we investigated the anzatz ( p.lO|) for the gold target and again found 



it to be a very good approximation at least for impact parameters smaller than the target radius. 

In order to be consistent with initial conditions ( |1.6|) we assume the following 6-dependence of 
AT^: 

N''{r±,x,xo;b) = (1 - e-''''^"'""'^^)^^)^ (2.12) 

with 

«:''(x,Xo,r^) = - ln(l - ,jND{r^,x,Xo)). (2.13) 

A^'^(r_L, X, xo) represents a solution of the equation ( |1.5| ) but with no dependence on the forth 
variable. The initial conditions for the function N^{r±, x, Xq) are set at 6 = and k^{xo, Xq, r±) = 
K{xo,r±). Since in the present paper we do not intend to compute cross sections, for which we 
would need to perform the b integration, the accuracy of the anzatz (p.l2|) will not be investigated 
here. 
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For each initial value of Xq the function N^{r±,x,XQ) is obtained after about ten iterations. The 
Fig. |I] shows the solutions as a function of the distance for various values of Xq and x. The 
amplitude for the elastic scattering iV^ |]16[ is plotted in the same graph. The obtained numerical 
inequality N'^ < < N is in perfect agreement with the physical expectations for the diffractive 
dissociation cross section to be larger than the elastic cross section. Another consistency check 
is the saturation of the function which is a consequence of the unitarity bound. In the black 
disk limit diffractive dissociation is a half of the total cross section. 



X = 10"^ 
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Figure 1: The function is plotted versus distance. The curves correspond to different values of xq: 



O - Xq = 10 ; b - Xq 



10 



-3 . 



c - xo = 10 ^. The solid line is N"^ . 



It is worth to investigate the dependence of the solutions obtained on the gap variable xq. To 
this goal we plot the function A^^ as a function of the gap Yq for various distances and at fixed 
y = 10 (Fig. 0). At short distances the solution depends strongly on xq though as we approach 
the saturation region this dependence dies out. 

It was stated in the Introduction that the function A^^ equals 2N — F, where both functions 
A^ and F are solutions of the same non-linear equation |T2|. Thus it is natural to compute 
2N — F solving the non-linear evolution equation |[Tl|, with appropriate initial conditions. A 
comparison with A^"^ from ( |1.5| ) would be an ultimate test for the correctness of the numerical 
procedures. Such test was successfully performed and we found an absolute agreement (relative 
error less than 1%) between both the computations. 



4 




Figure 2: The function N is plotted versus the gap Yq at fixed rapidity Y = 10. 



3 Saturation Scale 



Determination of the diffractive saturation scale Q^{x,xo) from the solution A^^ is a subject of 
this section. Unfortunately, no exact mathematical definition of the saturation scale is known so 
far. In the Refs. |]16|, several definitions of the saturation scale Qs{x) were proposed which 
related saturation scale to the shape of the function N. It is important to stress that it is not 
clear a priori whether should coincide with Qs or not. We will proceed here in the same spirit 



as Ref. |25|. Namely, we propose several definitions of the saturation scale while the variety of 
the obtained results will indicate the uncertainty in the definitions. 

For a step like function it is natural to define the saturation scale as position where it reaches 
half of the maximum: 



• Definition (a): 

iV^(i?f, x, xo) = 1/2, gf = 2/i?f. (3.14) 

The equality between the saturation radius and the saturation scale is motivated by 
the double logarithmic approximation. Though this approximation is formally not justified, we 
still believe it to make reliable estimates provided is large enough. The definition ( |3.14| ) is 
analogous to the one proposed in Ref. [jl6l N{2/Qs,x) = 1/2. If we recall that = N"^ at 
X = xq and postulate Q'^{xq,xo) = Qs{xo) then consistency requires 

• Definition (b): 

iV^(2/gf,a;,xo) = 1/4. (3.15) 

An alternative definition of the saturation scale could be one motivated by the Glauber- Mueller 
formula: 



• Definition (c): 

«^(2/Qf,a:,xo) = 1/2. (3.16) 

The saturation scales deduced through the above definitions are depicted in Fig. ^. For given 
Xo the observed hierarchy between the saturation scales obtained is an obvious consequence of 
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the definitions (|3A^ ), (|3J[5|) , (|3l6D and the shape of the function A^^ (Fig. |l|). Note that the 
saturation scale is almost xq independent. 

It is important to learn about x-dependence of the saturation scale. To this goal, we assume the 
following parameterization: 

Qf(a;,xo) = Qfo 4 ■ (3-17) 

In fact, the parameterization ( |3.17D is a good approximation for the values of the saturation scales 
obtained with 

A = 0.385 ± 0.015; and /? = 0.045 ± 0.025 . 

Within the errors these powers coincide for all the saturation scale definitions (|3.14| ), (|3.15| ), 
(|3.16| ). The small value for the power /? is a numeric indication of the very weak xo-dependence 
of the saturation scale. Its large relative error results on one hand from numerical limitations 
and on the other hand, this error signals for more complicated xo-dependence than it is given in 

& 

It is important to stress that the obtained power A coincides with the corresponding power of the 
saturation scale Qg [^ . 



4 Scaling phenomena 



In the Ref. p5| the function N was shown to display the scaling phenomena. We present here a 
similar analysis for the function N^. In the saturation region the scaling implies the amplitude 
to be a function of only one variable r = (r^ ■ (x, Xq))^: 



A^^(r^,x,Xo) 



(4.18) 



Let us define the following derivative functions assuming the scaling behavior ( |4.18|) : 



N^{r±,x,xo) 



dN^ dN^ (91n(g 



)D\2 



dY 



dr 



T ■ 



dlnx 



(4.19) 
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N^{r±,x,xo) = r[ 



ON 



D 



dN 



D 



dr ' 



dN 



D 



dYo 



dN^ d\n{Q 
dr 



D\2 



(9 In xq 



(4.20) 
(4.21) 



If the scaling behavior ([4.18|) takes place indeed, then both the ratios N^/N^ and ^/N^ are 
r± independent functions. Let us first consider scaling with respect to x. Fig. H presents the 
derivatives Nl^ and Nl^ as functions of the distance r±_ at fixed xq = 10~^. Both functions 
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2 4 6 8 10 
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Figure 4: The derivative functions (dashed line) and Ny (solid line) as functions of the distance 
at fixed xq = 10~^. 

and have extrema placed at the same distance depending on x. This is a consequence of 
the scaling behavior ( |4.18| ) and equations ( |4.19| ) and ( |4.20| ). The extrema occur at certain Tmax^ 
such that N'^'iTmax) = -TmaxN^"{Tmax)- In Fig. |, T^ax is approached by varying r± at fixed x. 
Alternatively it can be reached by varying x at fixed r± (Fig. |^). 




7 8 9 10 11 12 13 

Y 

Figure 5: The derivative function as function of the rapidity Y at xq = 10^^. 
Consider now the ratio function i?f : 

dHQi 



Ra{r±,x,Xo) 



y 



dlnx 



(4.22) 



If the scaling phenomenon takes place the function is expected to be r± independent. We 
study the scaling within the distance interval 0.04 GeV~^ < t± < 10 GeV~^ that corresponds to 



7 



0.25 GeV^ < < 2.5 x lO^GeV^. Fig. |^ presents the results on the scahng. The three hnes 
correspond to functions and divided by their minimal values within the interval, and the 
function i?f multiplied by the factor 40 to be seen on the scale. 




r^(GeV"') r^(GeV"') r^(GeV'') r^(GeV"') 



Figure 6: The scaling as a function of the distance at fixed xq = 10 ^. The positive curves are 
NP /N^^^^ (dashed line) and N^/N^^,-^ (solid line). The dotted line is 40 x . 

The function i?^ is clearly observed to be a very slowly varying function of r±_ for all values of x 
and r±. Though at fixed x the function cannot be claimed to be exact constant, its variations 
with r_L are very much suppressed comparing to the variations of the functions and N^^ . For 
example, at x = 10~^ within the given interval the function i?^ changes by maximum 20%, while 
within the very same interval both functions and change in several times. Then the 
relative fluctuation is much less than 10%, which confirms the scaling. The phenomenon holds 
with a few percent accuracy and it improves at smaller x ~ 10~^ and in the deep saturation 
region. However to observe this scaling behavior in these regions is numerically more problematic 
since both derivatives and tend to zero. 

The above analysis was performed for the fixed value xq = 10~^. Within the errors the function 
i?^ ~ —0.75 ± 0.08, constant independent on both rj_ and x. Moreover, if we repeat the same 
program but for different values of xq we discover quite similar scaling phenomena with being 
numerically independent on xq as well. This observation implies 

Qf(x,Xo) = Qfoixo)x~^; X = 0.37 ± 0.04. (4.23) 

Note that the value obtained for A is in agreement with the one determined in the previous section. 

Let us now study the scaling behavior with respect to the variable Xq. To this goal we investigate 
the function 3? which is related to the xo-dependence of the saturation scale Qfo- Assuming 
Qfg ~ Xq we predict 3? to have a maximum at r = r^m. Fig. |^ displays the function 3? as a 
function of the distance at fixed x = 4.54- 10"^ {Y = 10). In complete agreement with the scaling 
assumption ([4.21|) the function 3? possesses maximum with respect to r± variations. The heights 
of the maxima are proportional to (3. Since t N^'{T)\r=T^^^ — 0.2, (3 can be estimated to be 
approximately 0.05 ± 0.02 which agrees with the value deduced earlier. 

We can learn more about the scaling if we consider the function 3ft as a function of Xq or the 
energy gap Yq. In Ref. |2^ a model was built in which the function 3ft had a maximum with 
respect to Yq variation at fixed Y. We know now that this maximum is a consequence of the 
scaling phenomena. The dependence of 3ft on F = 10 is plotted in Fig. ^. 
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Figure 7: The function 5ft versus distance at fixed Y = 10. 
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Figure 8: The function 3? as a function of the ratio Yq/Y at fixed Y = 10. 



No maxima is observed on the plots of Fig. |^. In fact this is a sign of the scahng violation so 
far avoided by the discussion. The scaling with respect to Xq is not exact at Yq ~ Y. Due to its 
smallness (3? oc (3) the function 3? is most sensitive to small deviations from the scaling behavior: 



iV^(rx,a;,Xo) = N^^aUngir) + 5iV^(r^, a;, Xo) (4.24) 

In the kinematic region of the investigation variations of the function 6N^ with respect to rj_ and 
X are small compared to variations of N^^^^g. In contrary, the derivative of ^A^'^ with respect to 
Yq is of the same order as derivative of N baling- This is the origin of the large errors of P and the 
Xq scaling violation at Xq — x. 

In order to complete the analysis we propose yet another definition of the saturation scale based 
on the above presented scaling analysis. It is natural to define the saturation radius at the position 
where t N'{t) has maximum, namely at Tmax'- 



Definition (d): 

( d{rN'{r)y 



0. (4.25) 



4=4/(Q?)' 



The saturation scale obtained from ( |4.25 ) is depicted in Fig. Note again the weak dependence 
on the value of xq. 
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Figure 9: The saturation scale deduced from ( ^.25 ). 

The different curves correspond to xq = 10~^ 
(the upper curve), 10~^ (middle curve) and 
xq = 10~^ (the lowest curve) 
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5 Conclusions 

The non-linear evolution equation ( |1.5| ) is solved numerically by the method of iterations. The 
solutions obtained are in agreement with the unitarity constraints: the diffraction dissociation is 
larger than just the elastic scattering but smaller or equal than half of the total. 

The diffractive saturation scale is estimated form the solutions of ( |1.5|) basing on four different 



definitions of the saturation scale. Though there exists a significant uncertainty in the absolute 
values of the scale its x-dependence is found to be the same as of Qs - saturation scale deduced 



from the non-linear equation for [|TT], |T2]. In fact this result is quite natural. The dependence 
of the saturation scale on x is entire property of the evolution equation and it should not depend 
on both initial conditions and saturation scale definition. The saturation scale is discovered 
to be almost independent on the minimal gap xq. 

The scaling phenomena with respect to all variables were studied in details. The scaling with 
respect to x is well established. It holds with a few percent accuracy in the whole kinematic region 
investigated. The discovered scaling should manifest itself in the experiments on diffraction, and 
hence it would be interesting to search for it in the F^fi^x, Q"^) / {Q^ S) experimental data {S stands 
for the target transverse area). 

The numerically observed small scaling violation shows up when we consider the scaling with 
respect to xo- This happens due to the weak sensitivity of the solutions to the variation of xq. 
As a result, the variations of the solutions with respect to xq are of the same order as the scaling 
violation. The scahng sets in at x -C xq but is violated at x ~ xq. 

The detailed analysis of the ratio between the total diffractive dissociation and the total DIS 
cross section will be presented in a separate publication |^|. Our preliminary computations show 



that this ratio happens to be independent on the central mass energy in agreement with the 
experimental data [Q. This independence can be traced back to the scaling property displayed 
by the amplitudes N and A^^ and to the fact that both saturation scales depend on x with the 
very same power A. 
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